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We calculate the cross section of the inverse beta decay process, i/^ + n ^ p + e, in a. magnetic 
field which is much smaller than m^/e. Using exact solutions of the Dirac equation in a constant 
^ ' magnetic field, we find that the cross section depends on the direction of the incident neutrino even 

I when the initial neutron is assumed to be at rest. We discuss the implication of this result for pulsar 

^\ . kicks. 
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1 Introduction 



It has been argued by various authors that asymmetric neutrino emission during the proto-neutron star 
I phase can explain the observed recoil velocities of pulsars. One possible mechanism for this, first suggested 

by Kusenko and Segre [|l|, relies on neutrino dispersion |^ and oscillation in a background magnetic 



Oh, field. This requires non-zero masses for neutrinos. The other mechanism, which does not crucially depend 



on neutrino mass, relies on asymmetric neutrino cross section in a proto-neutron star. This can occur 
if the magnetic field in the proto-neutron star is asymmetric, as discussed by Bisnovatyi-Kogan His 
calculations were later checked and modified by Roulet who calculated neutrino opacity in strong 
magnetic fields and concluded that sizeable modifications from the zero-field case can occur. 



> 

X. 

^ ^ In this paper, we calculate the cross section of the inverse beta-decay process in a magnetic field. The 

process infiuences the neutrino opacity and hence the neutrino surface from which the neutrinos can escape 
freely from the stellar interior. The distorted 'neutrino sphere' produces a net momentum difference of 
the outgoing neutrinos which are emitted from the two hemispheres. Considerable work has been done 
on the URCA processes which have neutrinos in their final states Q . An angular dependence obtained 
1^ in the differential cross section of these reations imply that the neutrinos are created asymmetrically 
with respect to the magnetic field direction. In order to calculate the net momentum thrust for neutrino 
emission from a proto-neutron star, one should find out the angular dependence of the neutrino opacity 
owing to the presence of the magnetic field. This is the subject of the present paper. 

The neutrinos are assumed to be strictly standard model neutrinos, without any mass and consequent 
properties. We show that asymmetric neutrino opacity can arise even in a uniform magnetic field. This is 
because the presence of the magnetic field breaks the isotropy of the background, and a careful calculation 
in this background reveals a dependence of the cross section on the incident neutrino direction with respect 
to the magnetic field. 

The paper is organized as follows. In Sec ^, we provide some background for the calculation. Most 
of this section is not new, but we provide it for the sake of completeness, as well as for setting up the 
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notation that would be used in the later sections. In Sec. we define the fermion field operator and show 
how it acts on the states in the presence of a magnetic field. Sec. ^ contains the calculation of the cross 
section, which shows the asymmetric term mentioned above. In Sec. ^, we discuss the asymmetric term 
and estimate the magnitude of pulsar kick that can result from this term. Sec. H contains our conclusions. 



2 Solutions of the Dirac equation in a uniform magnetic field 

The Dirac equation in presence of a magnetic field is given by 

= [a ■ [p ^ eA) + Pm]iP , (2.1) 

where a and P are the Dirac matrices, and A is the vector potential. We will work with a constant 
magnetic field B along the z-direction. The vector potential can be chosen in many equivalent ways. We 
take 

Ao^Ay^A.^Q, A, = -yB . (2.2) 
For stationary states, we can write 

^ = e-^*(^), (2.3) 
where 4> and x £^re 2-component objects. Using the Pauli-Dirac representation of the Dirac matrices, we 



can then write Eq. (2.1) as 



{E~m)4> = (T-{p-eA)x, (2.4) 
{E + m)x = (T ■ {p ~ eA)(j) . (2.5) 



Eliminating x, we obtain 



{E^-m^)(t) ^ [a ■ {p - eA)f (t> 

= [p2 + eVs2 + eB(2yp,-a,)] 0, (2.6) 



where in the last step we have used our choice of A from Eq. ( |2.2| ). It is to be understood throughout 
that when we should write the spatial component of any vector with a lettered subscript, it should imply 
the corresponding contravariant component. Writing now p — — iV and noticing that the co-ordinates x 
and z do not appear in the equation except through the second derivative, we can write the solutions as 

^ = e^v-x^f{y), (2.7) 

where f{y) is a 2-component matrix which depends only on the y-coordinate, and the eigenvalue of the 
operator px , which we also denote by the same symbol. We have also introduced the notation X for the 
spatial co-ordinates (in order to distinguish it from x, which is one of the components of X), and for 
the vector X with its y-component set equal to zero. In other words, p ■ X^ = pxX -f PzZ. 

There will be two independent solutions for /(y), which can be taken, without any loss of generality, 
to be the eigenstates of tr^ with eigenvalues s = ±1. This means that we choose the two independent 
solutions in the form 



Putting these back into Eq. ( 21 ) and using the dimensionless variable 



2 



= 0, 



we find that Fg satisfies the differential equation 

£ + 

where 

E'^ — rn? ~ p1 + eBs 

as = — . 

eB 

This is a special form of Hermite's equation, and the solutions exist provided 
0, 1, 2, • • •. This provides the energy eigenvalues 



and the solutions are 



E"^ = +pI + eB{2n + 1 - s) , 



(2.10) 
(2.11) 

2n' + 1 for n' = 
(2.12) 



where Hn are Hermite polynomials of order n, and Nn are normalizations which we take to be 

/ \ 1/2 

^" " n! 2nA 



(2.13) 



(2.14) 



We stress that the choice of normalization can be arbitrarily made, as will be clear later. With our choice, 
the functions /„ satisfy the completeness relation 



^/„(0/„(eO ^V^SiC- e.) = S{y - y.) 



(2.15) 



where £* is obtained by replacing y by y^, in Eq. (2.9). 

Let us now classify the solutions by the energy eigenvalues 



Ei 



rn + Pz + 2neB , 



(2.16) 



which is the relativistic form of Landau energy levels. The solutions are two fold degenerate in general: 
for s ^ 1, n = n' and s = — 1, n = n' — 1. In case of n = 0, only the first solution is available. The 
solutions can have positive or negative energies. We will denote the positive square root of the right side 
by En- Representing the solution corresponding to this n-th Landau level by a superscript n, we can 
then write for the positive energy solutions. 









/„-i(0 



(2.17) 



For n = 0, the solution /_ does not exist. We will consistently incorporate this fact by defining 

/_i(2/) = 0, (2.18) 
in addition to the definition of /„ in Eq. ( ^T3| ) for non-negative integers 



The solutions in Eq. ( 2.17 ) determine the upper components of the spinor solutions through Eq. (2.7). 
The lower components, denoted by x earlier, can be solved using Eq. (2.5), and finally the positive energy 
solutions of the Dirac equation can be written as 



(2.19) 
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where denotes the space-time coordinate. And Us are given by 







Pz 



En + m 

Mr,. 



( 



where we have introduced the shorthand 







Pz 



(2.20) 



M„ = V2neB . 



(2.21) 



A similar procedure can be adopted for negative energy spinors which have energy eigenvalues E = 
'En- In this case, it is easier to start with the two lower components first and then find the upper 



components from Eq. (2.4). The solutions are 
where 



(2.22) 



Pz 



y+iy,n,p^) 



En + m 




y-{y,n,p^) 



Mn 

En + m 

p 



En + m^"-i(^) 


In-liO 



(2.23) 



with 



(2.24) 



For future use, we note down a few identities involving the spinors which can be obtained by direct 
substitutions of the solutions obtained above. The spin sum for the [/-spinors is 



Pu{y,y*,n,p^) = ^Us{y,n,p^)Us{y*,n,p^) 



2{En + m) 



{m(l + fT,)+;5[| -^||75}/n(e)^n(f*) 
+ {m(l - (T,) +^1,75} /n-l(e)/n-l(e.) 

-hM„(7i + ij2)IniOln-li^*) + Mn(7l ^ l72)/n-l (0^« (^*) 



(2.25) 



where we have introduced the following notations for any object a carrying a Lorentz index: 



11 ~ 
11 ~ 



(ao,0, Ojfl^) 
(a^,0,0, ao) 



(2.26) 

For the sake of completeness, we give the spin sum for the l^-spinors as well, although it will not be useful 
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for the rest of the present paper. 

Pv{y,y*,n,j\) = ^Vs{y,n,p^)Vs{y,n,p^) 



2{E„ + m) 



{-m(l + a,)+;5|| -^l|75}/„(0/n(C.) (2.27) 

M„(7i + Ij2)lni0ln-lit) - Mnill ~ ^2) In-liO In(I*) 



3 The fermion field operator 

Since we have found the solutions to the Dirac equation, we can now use them to construct the fermion 
field operator in the second quantized version. For this, we write 



s=± n=0 



dpx dpz 
D 



Mn,p^)e-'P-''^Us{y,n,p^) + f}{n,p^)e'P-''^Vsiy,n,p^) 



(3.1) 



Here, fsi^jPy) is the annihilation operator for the fermion, and /j(n,y;^) is the creation operator for 
the antifermion in the n-th Landau level with given values of px and pz ■ The creation and annihilation 
operators satisfy the anticommutation relations 



(3.2) 



and a similar one with the operators / and P, all other anticommutators being zero. The quantity D 



appearing in Eq. (3.1) depends on the normalization of the spinor solutions, and this is why Nn could 
have been chosen arbitrarily. It can be determined from the fermion field anticommutation relation, 
which is 

[i:{X),ijHX.)]^=S'{X^X^) (3.3) 

for — X'^. Plugging in the expression given in Eq. to the left side of this equation and using the 
anticommutation relations of Eq. ( |3.2| ) , we obtain 

dpx dpz 



[^{x),^\x.)]^ = j2j2 J 



^2 [e-^P-^---*h-^P^(^-^-^Us{y,n,p^)Ut{y.,n,p^) 

+e^P^i--^)e^P^(^-^^)V,iy,n,pjV,^{y,,n,p^)) . (3.4) 

Changing the signs of the dummy integration variables px and pz in the second term, we can rewrite it 
as 

[^{XUHX.)]^ = J2I1 I e-''^(— *)e-^^(--)(L/.(y,n,p^)C/t(y.,n,p^) 

+Vs{y,n,-p^)V,^{y^,n,-p^yj . (3.5) 



Using now the solutions for the U and the V spinors from Eqs. ( |2.2C ) and ( |2.23 ), it is straight forward 
to verify that 

E {Usiy,n,p^)Ul{y^,n,p^) + Vs{y,n,-p^)VXy*,n, -j^)^ 

S 

(En + my J L 



where 'diag' indicates a diagonal matrix with the specified entries, and ^ and involve the same value 
of px- At this stage, we can perform the sum over n in Eq. using the completeness relation of Eq. 
(2.15), which gives the (5-function of the y-coordinate that should appear in the anticommutator. Finally, 
performing the integrations over px and Pz , we can recover the (5- functions for the other two coordinates 
as well, provided 



2En 



1 



1 



£;„ +m 1)2 (27r)2 ' ^^-''^ 

using the expression for the energy eigenvalues from Eq. ( ^.16 ) to rewrite the prefactor appearing on the 
right side of Eq. (3.6). Putting the solution for D, we can rewrite Eq. (3.1) as 



dpx dpz / E, 



- EE/^ V .... 

X [/. (n, p^)e-^''-^^ {y, n,p^) + f} (n, p^)e'f Vs (y, n, ; 



The one-fermion states are defined as 



(3.8) 



(3.9) 



The normalization constant C is determined by the condition that the one-particle states should be 
orthonormal. For this, we need to define the theory in a finite but large region in the x and z-directions, 
of lengths Lx and Lz respectively. This gives 



2tt 



\/ LxLz 



Then 



i:u{X) \n,p,) ^ e-'^-^^£/.(y,n,H) |0) , 



where ipu denotes the term in Eq. (3.8) that contains the ?7-spinors. Similarly, 



(3.10) 



(3.11) 



(3.12) 



4 Inverse beta-decay 

In this section, we calculate the cross section for the inverse beta-decay process + n ^ p + in a 
background magnetic field. The magnitude of the field is assumed to be much smaller than m^/e or 
rrip/e, where to„ and rup are the masses of the neutron and the proton. For this reason, we can ignore 
the magnetic field effects on the proton and neutron spinors. The electron spinors, on the other hand, 
are the ones appropriate for the Landau levels. Thus, we can write the process as 



iy,ik)+niP)^p{P') + e{i^,n'). 



(4.1) 



4.1 The ^-matrix element 

The interaction Lagrangian for this process is 



•Cint — ^J^Gf- 



(4.2) 
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where L = |(1 — 75) and ~ Gp cos 9c, 0^ being the Cabibbo angle. In first order perturbation, the 
S'- matrix element between the final and the initial states of the process in Eq. (4.1) is therefore given by 

Sf, ^ V2Gpl d^X (e{i(,n') | ^^(,37''^^'^) | ^e(fe)) {p{P') |?(p)7p(.9v + 5A75)V'(n) | "(^)) •(4-3) 

For the hadronic part, we should use the usual solutions of the Dirac field which are normalized within 
a box of volume V, and this gives 



P{P') i^(p)l^ii.9v + gAlb)ik(n) n{P) 



J(P'-P)-X 



= [u{p)iP'h^igV + 9Al5)uin)iP)] ,(4.4) 



using the notations £ = Pq and £' = Pq. For the leptonic part, we need to take into account the magnetic 
spinors for the electron. Using Eq. (3.12), we obtain 



[U^„){y,n\j^)j^Lu^,^){k)] , (4.5) 



Putting these back into Eq. (4.3) and performing the integrations over all co-ordinates except y, we obtain 



En' + m 



2u!V 2£V 28'V2En,L^L^ 



1/2 



M 



fi ■ 



(4.6) 



Here, implies, in accordance with the notation introduced earlier, the (5-function for all space-time 
co-ordinates except y. Contrary to the field-free case, we do not get 4-momentum conservation because 
the y-component of momentum is not a good quantum number in this problem. The quantity M fi is 
the Feynman amplitude, given by 



Mf^ = ^2Gp^u^p)iP')-f^'{gv + gAl5)uin)iP) 
using the shorthand 

% = Py ^ ^ y ^ Py 



dy e 



iqyy 



(4.7) 



(4.8) 



The transition rate is given by |5/ip/T, where the quantization is done in a large time T. From Eq. 



(4.6), using the usual rules like 



Sis + LO - £ - En') 



T 

— S{£ + uj - £ - En 
2tt 



5{P^ + K~P',~p'J\ = ^5{P^ + K-P',-p',), 



SiP. + k,-P^-p'J 



2ti 



S{P, + k,-P^-p'J, 



we obtain 



\SMVT = -{27T)%'{P + k-P'-p 



, En' + m 



M 



16' ' " 'V^uj££'En' 

4.2 The scattering cross section 

Using unit flux \/V for the incident particle as usual, we can write the differential cross section as 



(4.9) 



(4.10) 



du^V 



dp, 



(4.11) 
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where dp, the differential phase space for final particles, is given in our case by 



(4.12) 



Therefore 



da ^ V 



T (27r)2 



(27r)3 

647r2 ^ ^ ^ ' LoS£'Er,> 



M 



'b^dp'^dp',_d^p'. 



(4.13) 



The square of the matrix element, averaged over the initial neutron spin and summed up over the final 
spins, is 



M 



(4.14) 



where H^'^ is the hadronic part, which is unaffected by the magnetic field because of our assumption 
stated earlier: 

H^^ = A{gl + g\){P>^P"' + P^P"^ ~ g^'' P ■ P') + A{gl - g\)mr,m^g^^ + ^gygAie^^'^Pp'y^P, . (4.15) 
The leptonic part on the other hand, is affected by the magnetic field. It is given by 



dy \ dy^ e''"*^"^*) Tr 



Pu{y*,y,n',j^)j^^j^L 



To perform the y integrations, we use the result 



which gives 



where 



V eB \VeB 



27r 1 
erS ynjn' + m) 



(4.16) 



(4.17) 



(4.18) 





4| 











+2M„52"^ 



Qy 



/eB 



In-l 



I 71-1 

qy 



% 

/eB 



/eB 



(4.19) 



Thus, 



M 



{gl + g\){P ■ KP' -k + P' ■ KP-k)-{gl^ g\)m„mpk ■ A 



+2gvgA{P ■ A P' ■ k- P' ■ A P ■ k) 



1 



27r 



eB {En' + m) 



(4.20) 



We now adopt the rest frame of the neutron and choose the axes such that the 3-momentum of the 
incoming neutrino is in the x-z plane. We will also assume that <C mp for the range of energies of 
interest to us. In that case, 



2t: mnrrirj 



^ eB En' + m 



idv + 35a)Aow + (<?y - gl)k,A, 



(4.21) 



The electron momentum p' enters the expression only through A", which has only the and 3 
components non-vanishing. The components of P' does not appear in this expression because of our 
assumption made earlier, except p'y which appears in A. Thus the integrations over fc^, p'^ and p'^ can be 
done easily in Eq. (4.13). The integrations overpf^ and p'^ get rid of the (5-fmictions over the corresponding 
components of momenta. As for fc^, we refer to Eq. (2_^). Since the center of the oscillator has to lie 
between —\Ly and \Ly, we conclude that —^LycB < p'^ < ^LycB. Thus the integration over k'^ gives 
a factor LycB. Putting back in Eq. (4.13) and using V — L^LyL^, we obtain 



da 



G} 5{Q + uj-En.) 



An 



LuEr, 



(.9v + 3gi)wAo + {glr ~ g\)k,K dP^dp', , 



(4.22) 



where Q is the neutron- proton mass difference, m„ — m^. 

We next perform the integration over Py. In the integrand, it occurs only as the argument of the 
functions In and In-i- The functions /„ are orthogonal in the sense that 



da /„(a)/„'(a) 



'eB 6r, 



(4.23) 



First of all, this shows that the term containing M„ in the expression for A" in Eq. ( 4.19 ) does not 
contribute to the cross section. The other two terms do contribute for n' > 0, and we obtain 



eBG'^ h ' 

dan'>o - + - E„,)\{gl + ?>g\) + [gl - g\)- 



dp'. 



Writing the argument of the remaining (5-function in terms of p'^ , we find that the zeros occur when 



(4.24) 



= p'± = ±^(Q + - m2 - Ml, . 
Thus, after integrating over p'^, we obtain the total cross section to be 



eBGl 



Cn'>0 



-(3y + 35i) 



M. 



/f2 



(4.25) 



(4.26) 



For n' — 0, however, the result is different, because in this case even the second term of Eq. (4.19) 
vanishes owing to the definition in Eq. ( ^.18| ). Thus, for this case, performing the integral over Py in Eq. 
([4.22D, we obtain 



dan 



eBGl 
47r 



5{Q + uj- Eo) 



{9v + '-^9a)-{9v-9a)- 



dp'z 



This gives 



eBGl 
2n 



i9v + S9l)-{gv-9l)- 



Q + uj 



The total cross section is then given as a sum over all possible values of n', i.e.. 



(4.27) 



(4.28) 



a = a„ 



eBGl V^^ 



n'=0 



27r 



n'=0 



9n'{gv + ^9a) - ^n',o{9v " 9a) — 



(4.29) 



where gn' is the degeneracy of the n'-th Landau level, 1 for n' = and 2 for all other values of n' . And 
"-max the maximum Landau level allowed at the given energy, given by the fact that the quantity under 
the square root sign in Eq. (4.26) must be non-negative, i.e., 

1 



mt 



2eB 



(4.30) 
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In Fig. |l], we have plotted the total cross section as a function of the magnetic field. The spikes in 
this plot appear at values of the magnetic field for which the denominator of Eq. (4.29) vanishes for some 
n' . For field values larger than this, that particular Landau level does not contribute to the cross section. 
The spikes in fact go all the way up to infinity, and their finite heights in the figure is an artifact of 
the finite step size taken in plotting it. In a real situation, where the initial neutrinos are not exactly 
monochromatic, these spikes are smeared out. 




0.5 



1.5 2 2.5 3 3.5 4 
log (eB/m^) 



Figure 1: Total cross section as a function of the magnetic field, normalized to the cross section in the field-free case. 
The initial neutrino energy is 10 MeV. The solid and the dashed lines are for the initial neutrino momentum parallel 
and antiparallel to the magnetic field. 



It is instructive to check that the results obtained above reduce to the known results for the field-free 
case. The contribution proportional to gy — g\ appears only in ctqi and vanishes in the limit _B — > 
owing to the overall factor of eB present in Eq. (4.28). Contributions proportional to gy + ig^ also have 
the factor eB with them, but in this case we also need to sum over infinitely many states. This gives 



Q + LO 



Q + oj 



^{Q + Lof - m2 - 2n'eB 2^{Q + wY -m? 



(4.31) 



where n'^^^-^ is the largest integer for which the quantity under the square root sign is non-negative. For 
_B — !■ 0, the last term vanishes, and we can identify n'^^^ as the integer for which the denominator of the 
summand vanishes. Thus we obtain 



'/3 I ..■2 



Wv + 3ffl) 



Q + oj 



^{Q + uo f -m? - 2n'eB 



(4.32) 



which is the correct result in the field-free case. 



5 The asymmetric cross-section 

The calculation of the cross section for inverse beta decay process has been performed earlier by Roulet 
H. He assumed that the matrix element remains unaffected by the magnetic field, only the modified 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

COS 9 



Figure 2: The cross section as a function of cos 6, where 9 is the angle the initial neutrino makes with the magnetic 
field. From bottom to top, the values of the maximum Landau level are 0,1,2 and 10. 



phase space integral makes the difference in the cross section. The resuhs he obtained is the same as the 
term proportional to gy + 3g\ that we obtained. 

It is the other term, not found by earlier calculations, that introduces the most important feature of 



the cross section obtained in Eq. (4.29), viz., its anisotropy. This comes only from the n' — contribution, 
which depends on the direction of the incoming neutrino momentum However, for n' ^ 0, there is a 
cancellation between the two possible states in a Landau level which washes out all angular dependence. 
In a real situation, then, the asymmetry will come only from the n' — state and its amount will depend 
on the relative contribution of this state to the total cross section. If the magnetic field is so high that 
only the n' = state can be obtained for the electron, the asymmetry will be large, about 18%. For 
smaller and smaller magnetic fields, the asymmetry decreases with new Landau levels contributing. This 
is shown in Fig. 

This fact can have far reaching consequences for neutrino emission from a proto- neutron star. It 
has been discussed in the literature that the presence of asymmetric magnetic fields can cause asym- 
metric neutrino emission from a proto-neutron star, thereby explaining the pulsar kicks. However, our 
calculations show that such an asymmetry in the magnetic field may not be necessary in producing the 
asymmetry in neutrino emission because of the fc^-dependent term in the cross section. Below, we make 
a rough estimate of the momentum asymmetry produced by this term. 

The typical size of a neutron star is about i? w 10 km. Thus, when a neutrino arrives at a density 
where its mean free path is about R, it escapes from the star. Therefore the condition for the neutrino 
to escape can be written as 

"nO-=-^, (5.1) 

where n„ is the neutron number density. We already observed that a is direction dependent. Therefore, 
the value of n„ on the "neutrino sphere" depends on the direction as well, and the surface is no longer a 
sphere. Different values of ri„ will correspond to different temperatures. Thus, neutrinos will be emitted 
with different momenta in different directions. This can result in a kick to the star. 



To estimate the magnitude of the kick, let us abbreviate Eq. (4.29) as 



cr = eBGl{a + bcos0), (5.2) 
where 9 is the angle between the magnetic field and the neutrino momentum. If we now consider the 
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directions 9^0 and 9 ^ n, the difference in neutron density on the corresponding points on the neutrino 
surface is given by 

2b 



Ann = 



eBGla^R 



(5.3) 



neglecting corrections of order b/a. 

The neutron gas in a typical proto-neutron star can be considered to be non-relativistic and degnerate. 
The number density of neutrons is thus given by | [lO| 



37r2 



1 + 



(5.4) 



where is the Fermi mometum, and we have neglected higer order terms in the temperature. This gives 

dUn mi ( T 



dT 



PF 



(5.5) 



So the temperature difference between the points on the neutrino surface in the 9 = and 9 = n directions 
is 

3 PF 26 



AT = 



ml T eBGla^R 



The momentum asymmetry can now be written as 

Ak _ 1 4AT 



(5.6) 



(5.7) 



6 T ' 

where we have assumed a black body radiation luminosity (cx T^) for the effective neutrino surface. 
The factor 1/6 comes in because the asymmetry pertains only to Ve, whereas 6 types of neutrinos and 
antineutrinos contribute to the energy emitted. This gives 

Ak A pp b 



k ml T2 cBGWR ' 



To find pp, we use the leading term in Eq. (5.4) and estimate n„ from the equation 

P(l - Ye) 

nn = , 

mp 

where Ye is the electron fraction and p is the mass density. Taking Ye = 1/10, we obtain 

pp = 24:pI{^ MeV . 



where pn is the mass density in units of 10 gem . Putting this back in Eq. (5.8), we obtain 



Afc 



1/3 



^ — 27 p^^ B^^ ^McV^2 ' 

where Bu = 5/(10^"* Gauss) and TmoV = T/ (1 MeV). For uo ^ m^ which is the relevant case, 

6=^^^ = 9.3x10-2, 
27r 

using gv = 1 and gA = 1.26. The value of a will depend on n^g^^- If o^^ly the n' = level contributes, we 
obtain 

(5.13) 



(5.8) 

(5.9) 
(5.10) 
(5.11) 

(5.12) 



a = = 9.2 X 10-1 . 

2tt 



This gives 



Afc 



1/3 r)-lrp-2 



MgV • 



(5.14) 



Obviously, with reasonable choices of p, B and T, it is possible to obtain a fractional momentum imbalance 
of the order of 1% which is necessary for explaining the pulsar kicks. 
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Figure 3: The value puVe above which no transition can take place at zero temperature due to Pauli blocking. 

6 Comments 

We have calculated the cross section of the inverse beta decay process in the background of a uniform 
magnetic field of arbitrary magnitude, using exact spinor solutions in a magnetic field. Our result shows 
that the cross section is asymmetric, viz., it depends on the direction of the incident neutrino. We have 
shown how this might explain the pulsar kicks. 

In the calculation, we have not taken the Pauli blocking effect into account. However, Roulet Q has 
pointed out that this effect is not very important in the range of densities in question. This can be seen 
from Fig. ||, where we have plotted the critical density beyond which the final electron is totally blocked 
at zero temperature because the Fermi energy is higher than Q + cj. At low but finite temperatures, the 
effect will be somewhat smeared. But it has been shown Q that in the relevant temperature range, the 
distortion of the Fermi distribution from the zero-temperature one is small. 

Acknowledgements : We thank Debades Bandopadhyay for many illuminating discussions and for 
pointing out some important papers in the field. We also thank A. Y. Parkhomenko for bringing Ref. Q 
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